Denaturation Transition of Stretched DNA 
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We generalize the Poland-Scheraga model to consider DNA denaturation in the presence of an 
external stretching force. We demonstrate the existence of a force-induced DNA denaturation 
transition and obtain the temperature-force phase diagram. The transition is determined by the loop 
exponent c for which we find the new value c = Av —1/2 such that the transition is second order with 
c = 1.85 < 2 in d = 3. We show that a finite stretching force F destabilizes DNA, corresponding to 
a lower melting temperature T{F), in agreement with single-molecule DNA stretching experiments. 
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Under physiological conditions the thermodynamically 
stable configuration of DNA is the Watson-Crick double 
helix. The constituent monomers of each helix, the nu- 
cleotides A, T, G, C, pair with those of the complemen- 
tary helix according to the key-lock principle, such that 
only the base-pairs (bps) AT and GC can form Upon 
heating or titration with acid or alkali of double-stranded 
DNA, regions of unbound bps proliferate along the DNA 
until full separation of the two DNA strands at the melt- 
ing temperature T^; depending on the relative content 
of AT bps, Tm ranges between some 60 to 110°C 01 . 
The classical Poland-Scheraga (PS) model views DNA 
as an alternating sequence of intact double-helical and 
denatured, single-stranded domains (bubbles or loops). 
Double-helical regions are dominated by the hydrogen 
bonding of bps as well as base stacking, bubbles by the 
entropy gain on disruption of bps The PS model is 
fundamental in biological physics and has been progres- 
sively refined to obtain a quantitative understanding of 
the DNA melting process [j, Q . DNA denaturation can 
also be induced mechanically, by longitudinal stretching 
of single DNA molecules by optical or magnetic tweezers 
or atomic force microscopes [a,0,@,0]- At the transition, 
the plot of stretching force F versus mean DNA extension 
L exhibits a plateau at 60-90 pN 1^, 11 1 . Force-induced 
destabilization of DNA has become a valuable tool, e.g., 
to probe the interaction of proteins that specifically bind 
to single-stranded DNA, at physiological melting temper- 
atures T„(F) well below T„(0) of free DNA |l3|. 





FIG. 1: Stretched DNA in the PS model with bound segments 
B and denatured loops Q. The DNA is attached between O 
and L and subject to the stretching force F in s-direction. 
Perfect matching in heterogeneous DNA requires both arches 
of a loop to have equal length £. 



DNA melting experiments. 

We treat the chain in the grand canonical ensemble in 
which the total number N of bps and the end-to-end vec- 
tor L fluctuate. The partition function in c? = 3 becomes 



Z{z,F) 



In this Letter we consider the force-assisted denatura- 
tion transition of double-stranded DNA in the framework 
of the PS model (Fig. [1]). The thermodynamic state of 
the DNA molecule now depends on both temperature T 
and stretching force F. We find a bounded region of 
bound states in the {T,F) plane (Fig. [2]). The shape 
of the transition line implies that finite stretching forces 
F indeed lower the melting temperature Tm{F), and 
the calculated force-extension relations F{L) exhibit a 
plateau over a certain range of DNA extension L (Fig. [3]). 
Both observations are in agreement with force-induced 



Ar=l 



dUZca„(iV,L)z^Vxp(/3FL 



(1) 



with (3 = \/{kBT). Zcan{N, L) is the canonical partition 
function of a chain of N bps with fixed end-to-end vector 
L and z is the fugacity. We assume that the force F acts 
in the positive x-direction, and is the a;-component of 
L (Fig. [T|) . If bound segments and bubbles are indepen- 
dent, Z factorizes: 



Z(z,F) 




Bfir > Bn, 



(2) 
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the last term equaling rigi?/(l — Bfl). The alternating 
sequence of bound segments and bubbles with weights B 
and f2 in Eq. Q is complemented by the weight of an 
open end unit at both ends of the chain. Note that only 
one strand of the end unit is bound to the, say, magnetic 
bead, while the other strand is moving freely. 

We model a bound segment with k ~ 1,2, . . . bps as a 
rigid rod of length ak where a = 0.34 nm is the length of a 
bound bp in B-DNA [13]. For simplicity we assume that 
the binding energy Eq < per bp is the same for all bps. 
The statistical weight of a segment with fixed number k 
and fixed orientation is then cj*'' with oj = exp(/?e) and 
e = —Eq > 0. Assuming that k fluctuates with fixed 
fugacity z, and rotates around one end while subject to 
the force F (Fig.[T]), the statistical weight of the segment 
for fixed z and F becomes 



B{z,u;,F) = / dnexp 

k=l ^ o 



{l3Fx) (3a) 



1 fl-Luze-y\ , , 

2y \i- — ujzey J 

Integration in Eq. pap is over the unit sphere with area 
47r, and x — ak cos 9 where 9 is the polar angle between 
segment and x-axis. At F = 0, B{z, lo, 0) — loz/{\ — u)z) 
as found previously for the denaturation transition of free 
DNA [13]. Note that B is only wcU-dcfincd for wze^ < 1; 
in what follows we assume z < e^^ /lu. 

Denatured loops are considered as closed random walks 
with 2£ monomers, corresponding to £ broken bps. This 
loop starts at O and visits the point r after £ monomers 
(Fig. [1]). The number of configurations of a loop is 



n{£,r) ^ Coi2£)pe{r) 



(4) 



under this constraint, where Co{2£) counts the con- 
figurations of a loop of length 2£ starting at O and 
Pi{r) is the probability that the loop visits r after £ 
monomers. For an ideal random walk in c? = 3, Co(2£) ~ 
^2e£-3/2 -j^j^g connectivity constant) and ^^(r) ~ 

7^"3 exp[-A(r/7^)2] where A > 0, r = ]r], and 7^ = M^^'^ 
is the scaling length of the walk. The amplitude b is 
proportional to the persistence length of the walk. Thus, 
ri(£, r) ^ s^^~^ exp[— A(r/7?.)^] where s = We assume 
that r moves freely and is subject to the force F in the 
positive x-direction. The weight of an ideal random loop 
for fixed £ and F is given by the Gaussian integral 

n{£, F)= I d^r n{£, r) e^^^ = As^^'^ exp{ay^£) (5) 



where A is an amplitude, c — 3/2, and a — 6^/(4Aa^). 
Finally, we sum il{£, F) over £ with weight z^ to obtain 
the statistical weight for an ideal random loop 

oc 

n{z,F) ^ A^u'^£-''^ AUc{u), szexp(ay^). (6) 



Lic(u) = X^fci "^^""^ the polylog function 14], con- 
verging for |u] < 1 for any c. For u — I three cases exist: 
(i) c < 1: Lic(l) diverges; (ii) 1 < c < 2: Lic(l) con- 
verges but Li'^(-u)|^_^ diverges; (iii) c > 2: Both Lic(l) 
and Li'^(w)|^ converge. The limit u = 1 corresponds 
to the value Zm{F) = exp{—ay^)/s of the fugacity; thus, 
il{z,F) is only well-defined for z < Zm{F) and diverges 
for z > Zjn (F) . The statistical weight fig of an end unit 
modeled as ideal random walk may be derived in a similar 
way and one obtains flei^, F) — v4eLio(w). 

For free DNA it was found that the nature of the denat- 
uration transition is determined by the analytic behavior 
of Lic(u) at u = 1: for c < 1 there is no phase transition 
in the thermodynamic sense; for 1 < c < 2 the transition 
is second order, and for c > 2 it is first order [i, 13]. One 
finds c = 3/2 < 2 if the loops are ideal random walks. 
Self-avoiding interactions within a loop modify this value 
to c = 3i/ = 1.76 with v = 0.588 in d = 3 [ll]. In both 
cases the transition is second order. Self-avoiding inter- 
actions between denatured loops and the rest of the chain 
were found to produce c = 2.12 > 2, driving the tran- 
sition to first order 1^ T^. These results suggest that 
the inclusion of self-avoiding interactions generally shifts 
the loop exponent c to larger values, possibly effecting a 
change of the transition from second to first order. 

To see how c changes when self-avoiding interactions 
within a loop are included for the case > 0, we obtain 
the weights and fie for a self- avoiding walk for ^ oo. 



Then, Eq. ^ holds with Co{2£) 



2tn-dv 



being the 



number of self-avoiding loops with 2t monomers. The 
probability density Pi[r) scales as ^^(r) — TZ^'^g{r/TZ) 
where TZ = bf is the scaling length of a self-avoiding 
walk and g{x) a scaling function. The function g{x) is 
not known for a self-avoiding loop. In what follows we 
assume g{x) ~ x'^ exp[— Ax''] for x ^ oo where A > 0, ^ 
is an exponent, and 5 = 1/(1 — i^) is determined by an 
argument by Fisher [l^. This form of g{x) is consistent 
with pi{y) for a Gaussian loop {v = 1/2, = 0) obtained 
above. For the related linear self-avoiding walk starting 
at O and ending at r after £ monomers, the above form of 
g{x) also holds and 4> can be expressed in terms of known 
exponents 15, [l^. For the present case of a self-avoiding 
loop 5 = 1/(1 — I') still holds but (j) is unknown. However, 
we will see that drops out from the result for ^{z,F) in 
the limit ^ — » oo at > 0. The integral in Eq. ([5]) is no 
longer Gaussian, but can be evaluated using the steepest 
descent method at k = fiFbl^ — > oo. It turns out that in 
this limit the integral is dominated by values r/^" oo. 
With the above behavior of g{x) at x ^ oo we find for a 
self-avoiding loop [cf. Eq. ^] 

n{e, F) = exp (ay^/''^) (7) 

for K oo with the new loop exponent in c? = 3, 



4iy- 1/2 = 1.85 



(8) 
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FIG. 2: Transition lines /m = Finale as function of i = ksT /e 
for Q = 1, s = 5 for denatured loops modeled as (a) ideal 
random walks and (b) self-avoiding walks (cf. Fig. [Sj. 



Thus, with self-avoiding interactions within a denatured 
loop and F > the transition remains second order, 
but moves closer to first order compared to free DNA 
(with c = = 1.76 obtained within the same ap- 
proach). The ampHtude A in Eq. ([7]) is proportional 
to the cooperativity parameter cto 1 quantifying the 
initiation of a loop in a previously intact double strand 
in the PS model |i, such that also A < 1. More- 
over, a = 0.6... 1.7 « 1 using a = 6^/(4Aa^) obtained 
for an ideal random walk where = 2LpXss /i] here, 
Xss = 0.6 nm is the length of a base in single-stranded 
DNA |lfl] and values for the persistence length Lp for 
single-stranded DNA were found to range between 0.7 nm 
d and 2nm [l8|. 

Finally, we sum r2(^, F) over t with weight to obtain 
the statistical weight for a self-avoiding loop [cf. Eq. ([5])], 



sz exp(a?/ 



(9) 



where 



1/(2^) 



FIG. 3: Force-extension curves / = Fa/e as function of I = 
{L^)/{a{N)) at fixed ti < t2 < ts < U (cf. Fig. Eb). Open 
circles mark second-order transitions (c = 1.85 < 2). Inset: 
f(l) for c = 2.5 > 2 where the transition is first order. 



fugacity is now given by Zm{F) — exp{—ay^^'^)/s. The 
weight rie(z,F) for an end unit obtains similarly, the 
result being Eq. ([9]) with c replaced byC = 3/2-|-z^ — 27 = 
-0.232, using 7 = 1.16. 

Phase diagram. We now obtain the transition line be- 
tween bound and denatured states in the (T, i^)-plane in 
the thermodynamic limit N ^ 00. For given fugacity z 
the average number of bps (open and closed) becomes 



(N) = d\nZ{z,oj,F)/d\nz, 



(10) 



where we explicitly include the argument lu from Eq. Q 
in the partition function ([2]). If is set one has to choose 
a fugacity z such that N — (N) ; in this case z becomes a 
function of oj, F, and A^. We denote the value of z in the 
limit A^ — > cxD by z*(u},F) = lim z(lu, F, N). Similar to 

the case F = 13], z*{uj,F) is the lowest value of z for 
which expression (fTO|) diverges. In the bound state the 
divergence turns out to occur when the denominator in 
(l — Bfl) vanishes [see text below Eq. ([2])], implying 
z*{uj, F) to satisfy 

B[z*,uj,F]fl[z*,F]^l, bound state. (11) 

Conversely, in the denatured state the divergence occurs 
because dz^e{z, F) diverges, which implies 



z*(uj, F) — Zm{F) , denatured state , 



(12) 



0.15 in c? = 3. The critical 



where z,n{F) is the critical fugacity obtained above 
(which is independent of lj). Thus, starting in a bound 
state in the (T, i*")-plane and approaching the transition 
line by varying T and F, the value z*{uj, F) is determined 
by Eq. (fTTj) and increases until it reaches the value Zm {F) 
from Eq. (H^. At this point the denaturation transi- 
tion occurs. In the denatured state z*{uj,F) is given by 
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Eq. ^2^. Right at the transition both Eqs. ((TTj) and (1121) 
hold simultaneously. Using Q[zm{F),F] — Ay^^Lic{l) 
by definition of Zj^iF) this impHes A(/3Fa)~^Lic(l) = 
1/ B[zjn{F),u!, F], relating F and to, or, equivalently, the 
reduced force / — Fa/e and temperature t = ksT je, for 
the transition line in the (i, /)-plane. 

The shape of the transition line fm{t) depends on A, 
a, and s. Fig. [2^ shows fm{t) for A = 1, a — 1, and s = 5 
for the case that denatured loops are ideal random walks 
[9 = Q,v = 1/2). The transition line for the more realis- 
tic value A ^ 1 is also shown (here A = 0.01). The line 
fm{t) separates a finite region of bound states from an 
infinite region of denatured states. The point (toi / = 0) 
with = tmif = 0) corresponds to the traditional melt- 
ing transition for free DNA {F — 0). The line fm{t) for 
A — I contains a region in which fm{t) decreases with 
t, such that increased stretching forces / lower the melt- 
ing temperature tmif), corresponding to force-induced 
destabilization of DNA [l3|. Interestingly, for A — 0.01, 
application of a small stretching force / first increases t.„i 
[l3. |20|. Moreover, frn{t) vanishes for both t to (as 
\t - iol^/^) and i ^ (as a^^/^i^/^). This means that 
for given < /o < / max, where /max is the maximum of 
fm{t), the chain does not only denature at a large im(/o) 



but also at a small i„(/o), as indicated in [2l|. This 
behavior can be traced back to a balance of the terms 
(/3Fa)^ and /3Fa in Zm{F) = exp{-ay'^)/s and Eq. (f3b| . 
respectively [13]. For {PFa^ <C l3Fa, i.e., ksT > Fa, 
the melting transition at ^^(/o) is mainly driven by the 
entropy gain on creation of fluctuating loops, similar as 
for free DNA. For fc^T <^ Fa the transition at tmifo) 
is due to the fact that B[zmiF),uj, F] decreases with 
y = [3 Fa = j jt in the denatured state, due to the rapid 
decay of Zm{F) [cf. Eq. ^\ Fig. [2)3 shows the line 
jmit) for self-avoiding loops with A = 1 and c = 1.85. 
Note that Eq. ([9]) reduces to the known result for a free 
self-avoiding loop (y = 0) only ii 9 = 0.15 is replaced by 
9 — 0; this is not a contradiction since Eq. ^ is based 
on the assumption that k = f3bF£'^ is large. To include 
in Fig. [21d the behavior of fm{t) for / = Fa/e — *■ we 
use Eq. ^ with 9 = 0.15 for y > 1 and 61 = for y < 1. 

Force- extension relations. In thermal denaturation of 
DNA one measures the fraction & of bound bps as func- 
tion of T. From the partition function ([T]) the average 
number (M) of bound bps is (M) = d\nZ/d\nuj and 
e = {M)/{N) with {N) from Eq. Conversely, 
stretching experiments on DNA reveal its response to an 
applied mechanical stress. The mean of the component 
of the DNA extension along F is (L^) = p-^d\nZ/dF. 
The average extension per bp in units of the bp-bp dis- 
tance a is / = {Lx) / {a{N)). For comparison with exper- 
iments and simulations we calculate I in the thermody- 
namic limit {N) oo. Consider the Gibbs-Duhem rela- 
tion for the thermodynamic potential lnZ(z,a;,F) [j^ j: 
Nd\nz + Md\nuj + (3LxdF = 0. If is fixed one obtains 
dlnz+Qdlnuj-'rldy — where z is a function of w, F, and 



N. FotN oo, d\nz*+Qd\nu;+ldy = 0, z*{u;,F) being 
the fugacity for (N) ^ oo as discussed above; Q{uj,F) 
and 1{lo,F) are the bound bp fraction and reduced DNA 
extension in the same limit. For constant y = (3Fa (or 
2/ = 0) one finds 9 = -91nz*(cj, F)/51nw (so that 9 
in the denatured state due to Eq. as expected). For 
constant w, corresponding to constant t = ksT/e, we 
find 1(10, F) = -{af3y^dliiz*{oj,F)/dF. Based on this 
result Fig.[3]shows force-extension relations f{l) at fixed 
^1 < ^2 < ^3 < t4 for the case that denatured loops are 
self-avoiding random walks (cf. Fig. ^jp). The two sets 
of curves correspond to expansions of /(/) for small / 
and close to the transition, respectively. The curves f{l) 
display flattened regions close to the transition, in qual- 
itative agreement with experimental force-extension re- 
lations for DNA. These regions become less pronounced 
as t increases and vanish for t ^ to = tm{F = 0). A 
force-extension relation for the case c > 2, for which the 
transition is flrst order, is also shown (here c — 2.5). In 
this case 1(f) jumps discontinuously from a value L to a 
larger value Z+ at the transition. 

We have shown that a longitudinal stretching force F 
results in a reduced denaturation temperature Tm{F), 
corresponding to force-induced destabilization of DNA. 
For the loop exponent in the presence of a finite _F > 
we found c = 4i/ — 1/2 = 1.85, so that the denaturation 
transition remains second order, but with an increased 
exponent. It would be interesting to study how the value 
of c is modified when self-avoiding interactions between 
a loop and the rest of the chain are included [3, [l3| • 
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FA9550-05-1-0472 (AH and MGO). 
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